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RESEARCH MEMORANDUM 


SUMMARY OF PITCH-DAMPING DERIVATIVES OF COMPLETE AIRPLANE 
AND MISSILE CONFIGURATIONS AS MEASURED IN FLIGHT 
AT TRANSONIC AND SUPERSONIC SPEEDS 
By Clarence L. Gillis and Rove Chapman, Jr. 

SUMMARY 


Longitudinal-damping data in the form of the pitching-moment damping 
derivatives ^C^ + Cm^ are presented and summarized from NACA flight 

tests of rocket-propelled models and full-scale airplanes. The experi- 
mental data are compared with calculations and a discussion of each con- 
figuration is given. Detailed conclusions are precluded hy the lack of 
systematic configuration changes and the relatively unknown effects of 
various factors. A general comparison of the results revealed the fol- 
lowing trends and general conclusions. The contribution to the pitch- 
damping derivatives of wings having 45° or less sweepback is erratic in 
the transonic region and may be either positive or negative, leading to 
possible dynamic Instability of tailless configurations using such wings. 

No satisfactory method of calculating the damping of such- wings In the 
transonic region Is available and a smoothly faired curve between sub- 
sonic and supersonic Calculated values will probably be unconservative. 
Configurations having triangular and swept wings with approximately 
60 ° leading-edge sweep and rounded airfoil sections exhibited less vari- 
ation in damping in the transonic region. Calculated values of the pitch- 
damping derivatives for configurations with horizontal tails are more con- 
servative and give better agreement with experimental results if the 
distance from the trailing edge of the wing mean aerodynamic chord (instead 
of from the center of gravity) to the tail Is used to calculate the down- 
wash lag effect. The damping caused by dovnwash effects arising from lift 
on a forward surface due to pitching velocity and rate of change of angle 
of attack may be appreciable if the forward surface is at a large distance 
from the center of gravity. 
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INTRODUCTION 


Successful accomplishment of aircraft flight at transonic and 
supersonic speeds has served to place increased emphasis on the dynamic 
behavior of aircraft. The type of behavior referred to here concerns 
motions of the aircraft as a whole as compared to motions of component 
parts usually referred to as flutter. Dynamic instability of the lat- 
eral and directional modes of motion is a problem that has occupied the 
designer* s attention for some time. A trend toward deterioration of 
the damping of the longitudinal motion of aircraft operating in the 
transonic region has also been experienced. This decrease in damping, 
as measured by cycles to damp, occurs because of the high altitude at 
which such flights are usually conducted and because of an actual 
decrease in the dimensionless aerodynamic damping- in-pitch derivatives 
at transonic speeds for some aircraft configurations (refs. 1 and 2). 

The- present paper summarizes the available experimental data 
obtained by the National Advisory Committee for Aeronautics on pitch- 
damping derivatives measured in flight at transonic and low supersonic 
speeds. Most of the data contained herein were obtained from rocket- 
propelled models of airplane or missile configurations flown at the 
Langley Pilotless Aircraft Research Station at Wallops Island, Va. 

Some data are also shown for several piloted -airplanes obtained from 
flights at Ames Aeronautical Laboratory and the NACA High-Speed Flight 
Research Station at Edwards, Calif. 

Some of the damping information presented herein has been published 
previously as parts of genera.l longitudinal-stability investigations and 
some Is, as yet, unpublished. In the present paper the damping data are„ 
compared with each other and with values calculated by the usual methods. 
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SYMBOIS 

lift coefficient, L/qS 

pitching-moment coefficient, fiLtehia g moment 


aspect -ratio 

constant In stability equation 
mass moment of inertia about Y-axis, slug-ft^ 
radius of gyration in pitch, j/ly/m. 
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lift 

Mach number 

wing area (including area enclosed within fuselage), sq ft 

horizontal-tail area (total included area), sq ft 

forward velocity, ft /sec 

weight of body, lb 

wing span in equation for c, ft 

damping constant in stability equations 

local chord, ft 

t/2 2 

c 2 dy, ft 

) 

acceleration due to gravity, ft/sec 2 

tail length, measured from center of gravity of configura- 
tion to center of pressure of tail, ft 

mass of body, W/g, slugs 

free-stream static pressure, lb/sq ft 

dynamic pressure when used in equations for Cl and Cm 
and equations (2) and ( 3 )> 7pM 2 , lb/sq ft 


mean aerodynamic chord, § 

b 


time, sec 

reference axis through center of gravity of configuration 
perpendicular to plane of symmetry 

angle of attack, radians 

specific heat ratio (l.4o) 

•control-surface deflection, deg 
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• 6 angle of downwash, radians 

6 angle of pitch, radians 

a rate of change of angle of attack, da/dt 

9 phase angle, radians 

cd frequency of oscillation, radians/sec 



Subscripts : ' 

f forward surface, based on its own area and chord 

r rear surface, based on its own area and chord 

t tail , _ 


T 


trimmed, or mean value 
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TEST AND ANALYSIS PROCEDURE 


Experimental data presented were obtained in free flight by the 
free-oscillation technique. In this test method, the aircraft is dis- 
turbed from a trimmed condition, usually by means of a rapid elevator 
deflection and the resulting short-period oscillation is recorded as 
the elevator is held fixed. The method of analysis of these oscilla- 
tions to obtain the static and dynamic stability derivatives is ade- 
quately covered in several references (such as ref. 3 ). For the 
present purpose, only that portion of the procedure dealing with the 
pitch-damping derivatives is of interest. 


Making the usual assumptions of constant velocity, level flight, 
and linear aerodynamic derivatives, a solution of the two-degrees-of- 
freedom equation of longitudinal motion of an aircraft can be obtained. 
For any appropriate quantity such as normal acceleration or angle of 
attack, the solution is of the form 


bt 


a = Ce cos (at + cp) + cap 


( 1 ) 


The constant b In equation (l) defines the damping of an oscillation 
and in terms of the aerodynamic derivatives is given by 


b = - 


Sq 

2mV 


1 / - \2 
Cr o, " 2( c “q + Cn &)\£y 


( 2 ) 


Equation (2) may be solved for the sum of the damping derivatives to 
give 


m, 


+ ^TTl* 


4ly 


q “a Sqc 


^( b + Hr c 0 


( 3 ) 


From the flight tests, therefore, the sum of the damping derivatives 
Cjnq and may be determined if the damping constant b is measured 

and the lift-curve slope is known. The damping factor b is gen- 

erally determined from the envelope of the curve defined by equation ( l) . 
This envelope is determined from the flight record of the appropriate 
measured quantity which, for rocket models, is angle of attack. The 
equation for the damping factor is 


log e Axg/^l 
t 2 - t-L 


(*) 
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where Axp and Axp are the amplitudes measured from the mean value 
of a at times tp and t2. For rocket-model tests, normal accelera- 
tion and angle of attack were measured to provide the lift-curve slopes. 
In some cases for the full-scale airplane tests reliable angle-of -attack 
information was not available from the flight test so the lift-curve 
slope was obtained from wind-tunnel tests and equations - (l) and (4) 
would be written in terms of the lift coefficient. The test and analy- 
sis method described does not permit separation of the derivatives Cmq 
and C^. This is not a severe limitation, however, because the damping 

is always proportional to the sum of the two derivatives regardless of 
the flight conditions or mass characteristics. 

In free-flight tests, accurate measurement of the damping deriva- 
tive is difficult for several reasons. First, as shown by 

equation (2), the total damping is composed of the damping-derivative 
term and the lift-curv^ slope term, the relative magnitudes of which 
depend on the radius of gyration in pitch. The inaccuracies in 
c mq + c m6> as obtaihed by solution of equation (2), are proportional to 

the relative contribution of the C^ term to the total, damping, b. 

As an example , the term contributed as much as two-thirds of the 

total damping in some of the rocket models, and in a full-scale test of 
an airplane (ref. l) the Cj^ term contributed about one-half of the 

total damping. Present design trends indicate that the proportion of 
damping contributed by the term on future airplanes will more 

nearly approach the proportion for the rocket -models . Other factors 
which may contribute to inaccuracies in measuring C^ + are non- 

linear aerodynamic derivatives, disturbances due to gusts, and any other 
effects on the oscillation peaks which define the envelope of the curve. 
Gusts and nonlinear .aerodynamic derivatives usually appear as apparent 
changes in the damping coefficients. 


CALCULATION METHOD 


The experimental damping derivatives presented herein are compared 
with calculated values from theoretical investigations wherever appli- 
cable, or experimental test data where such are available. It has fre- 
quently been assumed that all the damping on conventional airplane 
configurations is caused by the tail. A damping moment derivative Cmq 
results from the additional angle of attack at the tail caused by the 
velocity of pitching about an axis through the center of gravity. A 
damping moment derivative Cm^ results from the lag In downwash at the 
tail surface due to the finite time required for the downwash discharged 
at the wing to reach the tail surface (ref. 4). The damping derivatives 
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arising from these concepts are given by the equation 

(<S + c ”4) t ' i 1 + fe)( 2 °s)IKt) 2 (5) 

where the downwash lag contribution is represented by the de/da factor 
in equation (5). A factor such as this assumes that the downwash lag 
term is characterized by downwash discharged at the center of gravity of 
the configuration. It appears that this downwash lag term might be more 
representative if it is assumed that the downwash is discharged at the 
trailing edge of the mean aerodynamic chord of the wing rather than at 
the center of gravity. Modification of formula (5) in accordance with 
this concept gives 

(s. + S* - i 1 + & t# °SiKv) 2 (6) 


where 2’ is the length from the trailing edge of the mean aerodynamic 
chord to the center of pressure of the tail. 

For airplanes in which the tail contributes the largest part of the 
damping, equation (6) is satisfactory for an approximate calculation. 

It obviously fails for a tailless airplane. For airplanes with swept 
wings the damping contributed by the wing may be of appreciable magnitude 
at all speeds, and in the transonic region the wing damping may be of 
primary importance for wings of any plan form. Theoretical studies 
(refs. 5 to 8) and experimental data (refs. 9 and 10) show that, at 
transonic and low supersonic speeds, the wing itself may be dynamically 
unstable. Calculations of the damping derivatives Bhould therefore 
include the wing even though the damping due to the tail may be the 
major factor. 

Additional increments in damping-moment coefficients arise because 
of the downwash on a rear lifting surface resulting from the lift on a 
forward surface produced by the pitching velocity and the rate of 

change of angle of attack The damping coefficients caused by 

t hi s effect may be calculated as follows: 
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where S and c are the quantities on which C m is based. A similar 
expression can he derived for ACi^. The sum of the two damping-moment 
coefficients is then: 


A ( c mq + 


GLa r £r £r £ f dsW 

Cr ctf 5 s 5 



(9) 


As will he shown later this increment in the damping derivatives is very 
small for conventional airplane configurations hut can be fairly large 
for canard configurations. If the forward lifting surface is at a large 


distance from the center of gravity, the quantity 
to a close approximation by the expression 


d«r 


cam be obtained 


de_ = £^f de 
d^£j c da 


The theoretical results used to calculate the pitch-damping deriva- 
tives for subsonic and supersonic speeds are given In references 5 to 8, 
11, and 12. References 5 and 11 consider only untapered wings but the 
results were assumed to apply to tapered wings also. Wherever possible, 
wind-tunnel measurements of tail effectiveness and downwash were used to 
calculate the damping caused by the tail. 


In calculating the derivative C ®a the simple downwash lag concept 
as given in equation (5) and modified in equation (6) was used at all 
speeds. Several refinements to these calculations have been considered 
in references 13 and l4 but for the present purpose theje refinements 
were neglected because of the uncertainties in the experimental results 
caused by large possible experimental errors. in some cases and the rela- 
tively unknown but important effects of oscillation amplitude. 


RESULTS AND DISCUSSION 


The results of the flight measurements and the calculations are 
shown in figures 1 to 25 for airplane and missile configurations. All 
information pertaining to one configuration is given in one figure; 
break lines on the drawings of the models indicate wedge and hexagonal 
airfoil sections. The details of the airfoil section for each model 
are given in table I. For consistency, damping coefficients in all fig- 
ures are based on the total area and mean aerodynamic chord of the wing 
and are given for the rate terms (q and &) in radians_per second. 

This method of presentation does not correspond to thaiTof the references 
from which the data were taken for some cases. 
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Scales to which the derivatives ( C ^ are plotted vary con- 

siderably between figures. It is important to note, as pointed out in 
reference 15 , that a comparison of the absolute values of damping coef- 
ficients for different configurations has little significance. When 
presented, the curve designated "calculated (tail)" is the value obtained 
from equation (6), using wind-tunnel values of C]^. and de/da. The 

calculated (tail) curve thus includes the downwash lag effect which is 
an indirect effect of the wing but acts on the tail. The direct damping 
of the wing, calculated from theoretical results for the conventional air- 
plane configurations was added to the calculated (tail) curve to give the 
curve designated "calculated (tail plus wing)." For the canard configu- 
rations, the calculated curve includes the damping increments given by 
equation (9) and the theoretical values of Cm^ given by references 5 
and 6. There are insufficient data available to permit estimation of the 
effects of airfoil shape, oscillation amplitude, or oscillation frequency. 
Reference 1 6 shows some effects of oscillation amplitude and frequency on 
the damping derivatives (c^ + Cn&) for a 45° delta wing at subsonic speeds 

For reference, table II gives the reduced-frequency range for all models, 
and it may be noted that the frequency ranges for most of the data dis- 
cussed herein are below the ranges investigated in reference 1 6 . 


Rocket-Propelled Models 

Model 1. - The experimental data for model 1 are contained in refer- 
ence 17 and, as shown in figure 1, exhibit a smooth variation through 
the transonic region. Theoretical values of C m ^ + and are 

also shown in figure 1 for the pertinent supersonic speeds. The positive 
value of at supersonic speeds decreases the Cm^ damping of the 

wing by approximately one-third. Although the experimental curve is con- 
siderably higher than the theoretical curve, the variations with Mach 
number are similar. The effect of the fuselage was investigated by the 
use of reference l8, which does not include the effect of the afterbody, 
and was found to be negligible. No data are available on the effects of 
the afterbody. 

Model 2.- Figure 2 shows the experimental data for two models of a 
tailless airplane configuration. The data for model 2 were taken from 
unpublished rocket-test results. No curve is faired since the data were 
at isolated Mach numbers and were obtained from sustained low-amplitude 
oscillations (Ax w +0.25°) at M = 0.91 and 1.24. The oscillation 
amplitude for the points at M = 1.34 was Ax k +1.0°. 

Model 2B (data previously unpublished) contained pulse rockets to 
disturb the model and, therefore, higher angles of attack were obtained 
in the test. The data obtained at higher angles of attack (Ax « ±l.Q°) 
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confirmed the data for model A at M = 1.34, which was for similar angle- 
of-attack variations. Model B also exhibited a continued oscillation of 
low amplitude (±0.25°) similar to that which gave the positive damping 
coefficients for model A. Analysis of these low-amplitude oscillations, 
assuming time to damp to one-half amplitude 'a§“ infinite,;, gave values 
of Cmg + Cm& as varying from 1 . 98 to 2.25. The values of damping 

obtained at higher amplitudes shown by the data points for'model B and 
low amplitudes shown by positive data points for model A. are in reality 
boundaries for the damping coefficients where the coefficients for the 
system depends with certainty on the amplitude of the oscillation and 
possibly on the value of the mean angle of attack. 

In order to calculate the theoretical damping it was necessary to 
approximate the actual wing plan form by one more amenable to calcula- 
tion. Three approximations were tried: a swept tapered wing obtained 

by extending the leading and trailing edges to the root and tip, a- tri- 
angular wing having the same aspect ratio as the actual wing, and a tri- 
angular wing having the same leading-edge sweep as the actual wing. The 
calculated curves are shown for all three approximations and best agree- 
ment was obtained between the calculated 52 . 5 ° triangular wing and the 
flight data. _ 

Model 3 .- Data are presented in figure 3 for two tailless models 
designed to have different wing flexibilities. These data are from a 
program instituted to investigate the effects of wing flexibility on 
longitudinal stability. The solid, line (model A) in figure 3 is experi- 
mental data for a 9 -percent-thick wing of laminated wood and metal con- 
struction having medium rigid characteristics and is labeled "least 
flexible wing. " The dashed line is for a wing having the same geometric 
shape and thickness but with thinner metal inserts, hence, model B had a 
more flexible wing and is labeled "flexible wing." The least flexible 
wing had a stiffness of approximately one-half that of a' solid aluminum 
wing of the same thickness and the flexible wing had a stiffness of 
approximately one-third that of a solid aluminum wing of the same thick- 
ness. These stiffness ratios were determined “by comparison of the influ- 
ence coefficients for the solid aluminum wing and the laminated wing. 

A calculation of the damping of the wing was made -at M = 0.7 

by the method of reference 11. The subsonic calculation"gave the value 
of Qniq as approximately - 1 6 , which is far too high and unconservative 

when compared with the experimental data. It is possible that the disa- 
greement between calculation and experiment at subsonic speeds may be 
due to a positive for the wing which is not accounted for in the 

computation. A computation of damping at M = 1.4l4 was made by using 
references 7 and 8 and good agreement with the experimental data was 
obtained. 
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Model 1)-.- Unpublished data for a 45° swept-wing tailless configu- 
ration are shown in figure 4 . The wing was a 6-percent-thick section 
constructed of solid dural. Qualitatively, the shape of the curves for 
damping of models 3A, 3B, and 4 show remarkable similarity for Mach num- 
bers from 0.90 to 1.15 with damping for both models exhibiting a charac- 
teristic drop in this region. The wing on model 4 was comparable in 
stiffness to the flexible wing of model 3. Data for model 4 were not 
obtained below a Mach number of 0 . 88 ; hence, it is not known if the damping 
decreased in a manner similar to that for model 3 . One calculated point 
for the wing alone is shown at a Mach number of 1.38* The calculated value 
is conservative but shows poor agreement with the experimental data. A 
calculation of body damping by the use of reference 19 showed the body con- 
tributiqn to be small at supersonic speeds because of the afterbody shape. 

Model 5 .- The data for model 5 (fig. 5 ) are from reference 20. 

This model was a wingless fuselage-tail configuration that has been used 
for tests of a number of wing plan forms for supersonic airplanes 
(models 6 to 12). The damping coefficients in figure 5 are based on the 
wing area and chord of model 6 for compar ison purposes. The theoretical 
values of C m g + agree fairly well with the experimental values, 

particularly at supersonic speeds, as does the calculated values of 
of the tail surface (curve labeled "calculated Cnig (tail)"). The agree- 
ment of the theoretical and experimental values 1 b somewhat fortuitous, 
however, because, as shown in figure 5> the theoretical values of 
are considerably higher than those calculated from measured values of 
tail effectiveness. The theoretical values of 0^, which are those for 
the isolated tail uninfluenced by downwash, are positive throughout the 
supersonic speed range considered and become rapidly larger aB the Mach 
number decreases below about 1.2. 

Model 6 .- The data for model 6 (fig. 6) are from reference 20. Two 
geometrically identical models were flown; model B having a steel wing 
(designated "rigid wing" in fig. 6) and model A, having an aluminum wing 
(designated "flexible wing" in fig. 6). At transonic Mach numbers rather 
large changes in the experimental damping data occurred as contrasted to 
the smooth variation with Mach number of wing-off data (fig. 5 ) and calcu- 
lated tail contribution. As indicated by the calculated curve including 
the direct wing effects, part of the loss in damping near M = 1.0 can 
be explained by a large positive value of for the wing. This 

large destabilizing effect of the wing is caused by the positive value 
of Cju^ being much larger than the negative value of C^. Another 

reason for the large damping changes near M = 1.0 can probably be 
found in unaccounted-for downwash changes in this region. In calcu- 
lating the Cm^ portion of the tail contribution to damping it was 

assumed, in the absence of appropriate data, that the downwash variation 
with Mach number was similar to the variation of lift-curve slope, which 
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did not show large changes (ref. 20). Since large changes in static 
stability for this configuration took place in the transonic region 
(ref. 20), it is probable that the dovnwash variation with Mach number 
was considerably different from that assumed. 

The increment in + Cn^ caused by the downwash resulting from 

the and on the wing was calculated for model 6 and is shown 

dn figure 6. This contribution to the damping is very small because the 
terms CLq and C]^ are very small and, at least for supersonic speeds, 
are of opposite sign. Since this condition should be approximately the 
same for all configurations having the forward surface near the center 
of gravity, this contribution to the da, raping was neglected for all other 
such configuratibns herein. 

Model 7 »- Model 7 was identical to model 6 except that the wing 
had an NACA 65(o6)-A-004.5 airfoil section instead of £he sharp-edge 
hexagonal section on model 6. Data for model 7 are shown in figure 7 
and were taken from- unpublished data. Model 7 , which had an aluminum 
wing, showed less damping than does model 6k and a larger disagreement _ 
between the low- lift experimental curve and the calculated curve for the 
tail. The reason for the disagreement in the three sets of experimental 
data has not been established. 

Model 8 .- Two experimental curves (fig. 8) were obtained for model 8 
(ref . ”21)": one at low lift coefficients and one at high lift coefficients 

The calculated tail damping shows a large decrease at Mach numbers near 
0.93, as did the test data, and this is caused by the large decrease in 
downwash at these Mach numbers (ref. 22). The direct effect of the wing 
was fairly small in the region where it could be calculated. The direct 
contribution of the wing (as computed) is destabilizing at the supersonic 
speeds of the test. It Is perhaps significant that when the effect of 
the wing is included, the variations with Mach number of the experimental 
and calculated curves are, in general, similar at supersonic speeds. 

Model 9 .- The damping data for two k^° swept-wing models are pre- 
sented in figure 9 and are not previously published. Model A had a 
steel wing and is denoted as "rigid wing" in the figure. Model B had 
an aluminum wing and the curve is labeled "flexible wing" In figure 9 . 
Calculated values of the damping shown for models A and B are based on 
experimental tail effectiveness. An additional curve labeled "calcu- 
lated (tail, theoretical)" is shown in the figure and is based on theo- 
retical tail effectiveness. The direct wing contribution to damping at 
pertinent supersonic Mach numbers was small and is not included in the 
calculated curves. It is noted that the difference in the experimental 
damping curves for models A and B is in the opposite direction, to that 
difference in the calculated curves for the two models, which is a 
center -of- gravity effect. A difference such as this might possibly be 
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explained by a higher de/da for the more rigid wing. Severe changes 
in downwash for the flexible wing might explain the high peak and rapid 
change that occurs between Mach numbers of 0.93 and 1.03 for the flexible 
wing. Peak values and severe changes such as these occur on a swept-wing 
airplane which is discussed in a subsequent section. 

Model 10 .- The data for a 60° swept wing model are shown in fig- 
ure 10 as obtained from reference 23. The direct contribution of the 
wing to the total damping is illustrated by the increment between the 
total damping curve and the tail-damping curve (includes downwash lag). 

Two additional curves are shown in the figure to illustrate different 
techniques of computation. The curve labeled "calculated (tail, theo- 
retical)" uses theoretical tail effectiveness plus a downwash lag term 
from estimated downwash. The peak value of damping which occurs at 
M = O .85 is not predicted by the calculations but is similar to a peak 
that occurs in the data of model 11. 

Model 11 .- The experimental data presented in figure 11 for model 11 
are reported in reference 2k. This modified triangular-wing model exhib- 
ited damping characteristics that showed a marked change between Mach 
numbers of 0.84 and ‘0.88 where the damping derivatives almost doubled in 
value in a manner similar to model 10. Calculated values for the tail 
and for the tail plus wing are also shown in the figure and show a smooth 
variation over the Mach number range. The difference in variation between 
the experimental and calculated curves and a comparison with the data for 
model 5 indicate that the rise in damping between M = 0.80 and M = 0.90 
can be attributed directly to the wing or to increased downwash from the 
wing. 


Model 12 .- The data presented in figure 12 for the 60° triangular- 
wing model are from reference 25. Two experimental curves are shown in 
the figure: one for the low- lift-coefficient range and the other for 

the high- lift-coefficient range. This difference in damping for high 
and low lift coefficients is primarily attributed to the nonlinear down- 
wash behind the triangular wing. The calculated (low lift, wing plus 
tail) curve shown is for de/da evaluated at zero lift and for experi- 
mental tail effectiveness. Calculated damping coefficients, shown in 
figure 12, for the wing-plus -tail (high lift) configuration indicate 
that the increased damping at high lift can be primarily attributed to 
the increased downwash. 

Model 13 .- Damping data for a supersonic airplane configuration are 
shown in figure 13. Experimental data are shown for two models; model A 
data were at high lift and model B data were at low lift. Curves for 
tail damping and total damping for model B are also shown in the figure. 
No computations were made for the high- lift model. The disagreement 
between the two sets of experimental data at subsonic speeds is believed 
to be due to the changes in downwash for high and low lift coefficients. 


* 
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Good agreement was obtained between the computed and experimental data 
at low lift. From the experimental curve and the fact that it fairs 
into the calculated (tail plus wing) curve, it appears that the direct 
wing contribution to damping approaches zero _at Mach numbers of approxi- 
mately 1.17 and 0.90 but generally does not give much destabilizing 
effect. It is noted that if the calculated (total damping) curve were 
extended to lower supersonic Mach numbers, the unstable contribution of 
wing would make the agreement between calculation and experiment 

worse. 

Model l4 . - Data for a 45° swept-wing airplane model from reference 26 
are presented in figure lb. Calculated curves are shown for the tail 
damping using both the modified and unmodified tail length for the down- 
wash lag term. Flight data vary unpredictably through the transonic 
region and a comparison with the calculated curve for the tail indicates 
that the wing contributes both positive and negative damping over the 
region. A single data point (unpublished) for a similar tailless con- 
figuration shows the wing to have a positive (unstable) damping coeffi- ~ 
cient at M = 0.79- The calculated curve using the modified tail-length 
concept gives better average agreement with the experimental data and is 
more conservative than the unmodified formula. No total damping calcu- __ 
lated curve is shown since the wing contribution was small at the higher 
Mach numbers where it could be computed and could not be computed over 
the Mach number range where data were obtained. 

Model 13 . - The experimental data shown in figure 15 for model 15 
were taken from reference 27. A curve of calculated tail damping shows 
fairly good agreement with the experimental data, both as to magnitude 
and general variation with Mach number. At subsonic speeds the calcu- 
lated wing damping was appreciable, being about one-third of that con- 
tributed by the tail. For the supersonic speed range below a Mach number 
of 1.28 the calculated wing damping coefficient is positive and the cal- 
culated damping of the complete configuration indicates a greater vari- 
ation with Mach number than is shown by the experimental data. 

Model 16 .- Reference 28 and figure 16 contain the experimental data 
for model l6. Data are presented for models A and B, identical models 
except for center-of-gravity location. The effect of the center-of- 
gravity location on the damping is small. In contrast to' most other 
swept-wing configurations included herein, the variation with Mach num- 
ber of the experimental damping coefficients Is smooth and the damping 
increases as the Mach number Increases. The significance of these 
phenomena and whether they are related to the inverse taper are not 
known. At subsonic speeds, the calculated wing damping is about 50 per- 
cent of that calculated for the tail. Because of the Inverse taper 
ratio on the wing, the methods of references 7 and 8 are not applicable 
to this model. 
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Model IT .- The damping data for model 17, shown in figure 17, are 
from unpublished test results. "This model had the same wing as the tail- 
less model 3 and had the same characteristic loss in damping near a Mach 
number of 1.0. Comparison of the experimental curve for model 17 with 
the calculated damping of the tail indicates that the wing contributed a 
positive ('unstable) damping coefficient near M = 1.0 as did model 3, 
although the quantitative values of wing damping for models 3 and 17 are 
not directly comparable because of the different center-of -gravity loca- 
tions. The agreement of the total calculated damping with the experi- 
mental values is very good at subsonic speeds. The wing damping could 
not be calculated over the supersonic speed range of the tests. 

Model l8 .~ The results shown in figure l8 for a 60° triangular -wing 
missile configuration are from reference 29. Damping data are presented 
for both the inline- and interdigitated-tail positions. As pointed out 
in reference 30 the marked difference in the characteristics for the 
inline and interdigitated tail can be attributed to the difference in 
downwash at the tail by considering the rolled-up vortex sheet. The 
experimental curve for the inline-tail model has a variation with Mach 
number that is, in general, similar to that for the modified- triangular- 
wing model (model ll) . A calculated curve for the tail is shown as well 
as a calculated curve for total dangling coefficient (both curves being 
for the inline-tail configuration) . Estimates of downwash and tail 
effectiveness based on experimental data were used and good average 
agreement with the flight data were obtained. A theoretical curve is 
shown and is generally nonconservative at the higher Mach numbers but 
shows good agreement with experimental data in the region near a Mach 
number of 1.1. 

tfodel 19 . - The damping data for a canard missile model shown in 
figure 19 are unpublished. The calculated damping is considerably less 
than the experimental damping, but the variation with Mach number of the 
two curves is remarkably similar. The increments in damping coefficients 
caused by the downwash resulting from the canard lift effects 

and were calculated and found to be appreciable as shown in fig- 

ure 19. This was true for all the canard configurations considered 
herein. The reason is that the value of Cl^ depends upon the distance 

of a lifting surface from the center of gravity (refs. 5 and 6) and, for 
the canard surface, becomes fairly large and of the same sign as Cj^, 

resulting in an appreciable effect on the damping. 

Models 20 and 21 . - The data for models 20 and 21 were taken' from 
references 31 and 32 and are shown in figure 20. Models 20 and 21 were 
geometrically identical except that the cylindrical portion of the body 
between the wing and the canard was 40 inches longer for model 21 than 
for model 20. All the experimental curves in figure 20 show similar 
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variations with Mach number near a Mach number of 1.0, 'which differ con- 
siderably from the variations indicated by the calculated curves. For 
model 21 neither the theoretical damping nor' that calculated from experi- 
mental data predicted the variation with Mach number of^the measured 
damping a,t supersonic speeds. The calculated' values for model 20 were 
considerably lower than the experimental values at supersonic speeds but 
showed a similar variation with Mach number. 

Model 22 .- The experimental values for model 22 (fig. 21) were 
taken from unpublished results. Models A and B were identical except 
for center-of-gravity location and the different lift-coefficient ranges 
covered by each model. Calculated values exhibit fair agreement with 
measured values only at supersonic speeds and low lift .coefficients. 

The difference between the calculated values .for models A and B was 
caused partly by the different center-of-gravity locations and partly 
by different values of downwash (shown by wind-tunnel data) over the two 
lift- coefficient ranges. Also shown in figure 21 are tbe values of 
damping contributed by the downwash resulting- from the canard lift 
effects and " t* 


Full-Scale Airplanes 

Airplane 1 .- Damping data from flight .tests of the jC-1 research air- 
plane are shown in figure 22. These data were obtained from reference 33 , 
and reference 22 was used to obtain the calculated curves shown on the 
figure. The calculated damping agrees well with the experimental data 
both as to magnitude and variation with Ifech number. The calculated 
total damping curve is somewhat lower and agrees better with the expert- _ 
mental data than the calculated curve given in_ reference 33- Thi6 is 
probably caused by the use herein of the modified tail length for com- 
puting the downwash iag effect (eq. (6)). Model 8, previously discussed, 
was equipped with the X-l wing. The low-lift ^experimental data for 
model 8 show a variation with Mach number very similar to that for the 
flight-test data of the X-l airplane. 

Airplane 2 .- The data for the X-4 airplane are shown in figure 23, 
as obtained from reference 2. Measured values are considerably lower 
than the calculated values and, for the highest altitude "test, became 
zero at a Mach number of about O. 85 . At a Mach number of 0.88, a sus- 
tained oscillation of small amplitude was obtained in the flight test 
which yielded a large positive value of + C^. Because this oscil- 

lation was of small amplitude, the value shown can be considered valid 
only for small- amplitude motions. The measured values showed an effect 
of altitude, which, as stated in ref erence 2, may dedicate that aero- 
elasticity has a large effect on the damping in this case. 
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Airplane 3 .- The results for the F-86 airplane in figure 2k were 
taken from reference 1. The sharp drop in the experimental curve at 
M = 0.92 was confirmed by a number of separate test points and corre- 
sponds to a similar irregular variation in the static stability 
parameter Cn^. Definite reason for these fairly sharp changes was not 

established in the preliminary investigation of reference 1 but they are 
similar to those for other swept-wing configurations (models 3, 4, 9 , 

.and 15 ) . 

The estimated damping contributed by the tail and the tail plus 
wing are shown in figure 2k. Tail contributions were calculated from 
lift-curve slopes and downwash values obtained from transonic bump tests 
of 35° swept wings. Although reference 34 contains tail effectiveness 
data (dCm/dit) these were not used in calculating the damping because 
the values measured were about 25 percent lower than would be obtained 
by using the lift-curve slopes of an isolated 35° swept-tail surface. 

This is probably primarily caused by the fact that the inboard portion 
of the adjustable stabilizer (approximately 10 percent of the span) does 
not move. The lift contributed by the entire surface, however, is 
effective in producing damping. The damping contributed by the wing was 
appreciable at subsonic speeds for this configuration and a comparison 
of the calculated and experimental data indicates that the damping due 
to the wing is approximately zero at Mach numbers near 1.0. 

Airplane 4 .- The damping data shown in figure 25 for the D-558-II air- 
plane were taken from reference 35. The calculated damping agrees fairly 
well with the experimental data at Mach numbers up to about 0.75, but as ■ 
the Mach number increases, the experimental data Bhow, first, a more rapid 
increase, and, then, a more rapid decrease than the calculated values. 

The more rapid decrease may be attributed, at least partially, to a proba- 
ble dynamically unstable contribution of the wing as shown by model 3 and 
airplane 2. 


GENERALIZATION OF RESULTS 


A generalization of results on damping derivatives is limited by 
the number and types of configurations and the continuity of the experi- 
mental flight data available for summary. Some general trends are noted 
for each type of configuration and a discussion is presented on the com- 
putation of damping. The general trends noted are limited for some of 
the cases by inconclusive evidence. 
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Tailless Configurations 

Four of the five tailless configurations (models 2,-3, and 4, and 
airplane 2) exhibited either very small negative values or positive 
values of Cm^ + C^ in the transonic region. The fifth configuration 

(model l), with a 60° triangular wing, showed a fairly smooth variation 
through the transonic region with no indication of a decrease in damping 
similar to that encountered on the swept -wing configurations. On model 2 
and airplane 2 the positive (unstable) damping coefficients were obtained 
from small- amplitude sustained oscillations and the. evidence indicates 
that the damping is nonlinear with amplitude, being less. at the smaller 
amplitudes. The data shown for model 1 were obtained from fairly large- 
amplitude oscillations but other models having the same wing on a slightly 
different fuselage have been flown with no pitch disturbance applied 
(unpublished) and no sustained oscillations were encountered. 

Examination and comparison of the data for models 3A, 3B, and 4 
show that the damping of swept wings is dependent on the. flexibility. 

Also, these data show that similar swept wings exhibit similar varia- 
tions of damping with Mach number, characterized by a rapid decrease in 

damping above M = 0.85 which would not appear to be a critical Mach 

number effect. This same decrease in damping is evident! for airplane 2. 


Wingless Configurations 

Data are shown for only one wingless configuration (model 5) and 
indicate little variation in the damping coefficients with Mach number. 
It is reasonable to expect that the damping contributed by any lifting 
surface having a large moment arm from the center of gravity should vary 
with Mach number in a. manner similar to its lift-curve slope if it is 
not operating in the flow field from a forward surface. 


Wing-Plus-Tail Configurations 

None of the configurations having tail surfaces exhibited any posi- 
tive values of C^q + Cm£ through the Mach number rangejtested, although 

several experienced large changes in the transonic region resulting in 
some fairly small negative damping coefficients. The contribution of an 
unswept wing to the damping derivatives of the total configuration is 
usually erratic and may be either positive or negative in the transonic 
region. In the majority of the cases the damping of the unswept-wing 
configurations was nonlinear with lift coefficient. 

In general the swept-wing configurations having 45° or less sweep 
also showed erratic changes in damping in the transonic region although 
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such changes were not as severe as those for the tailless swept-wing 
configurations because of the smaller contribution of the wing to the 
total damping. The one 60° swept-wing configuration discussed had a 
fairly smooth variation with Mach number near M = 1.0. The triangular- 
wing configurations generally had a somewhat smoother variation with Mach 
number than those with unswept or swept wings. Modifying or thickening 
the airfoil sections of the triangular wing (models 11 and 18 ) seems to 
result in somewhat less smooth variations in damping with Mach number. 

The results for model 17 show that the orientation in roll of the 
cruciform wings and tails with respect to each other had a significant 
effect on the damping. It has previously been shown in reference 36 
that the tail bending of a configuration similar to model 18 appreciably 
affected the damping in roll. It is probable that tail flexibility 
effects also influence the pitch damping. 


Canard Configurations 

It appears that for the canard configurations shown, the experimental 
data are less amenable to accurate prediction, probably because of large 
interference effects and in the case of model 22 by nonlinearities in all 
the aerodynamic characteristics. Practically no experimental data are 
available at transonic speeds on the downwash and interference effects 
applicable to canard configurations. 

The contribution to the damping coefficients of the downwash caused 
by the canard lift increments and was appreciable for all the 

canard configurations considered and it appears that this effect should 
not be neglected for any configuration in which the forward lifting sur- 
face is at a large distance (either forward or rearward) of the center 
of gravity. 


Agreement of Calculations apd Experiment 

It has been generally shown that a computation of the damping con- 
tributed by the tail can be made and It will agree reasonably well with 
the experimental damping. If it is assumed that the downwash lag effect 
of the wing on the tail is that effect caused by the change in downwash 
in the finite interval of time required for air to travel from the 
trailing edge of the mean aerodynamic chord to the tail rather than from 
the center of gravity to the tail, then the damping due to the tail will 
exhibit better agreement with experimental data and will be a more con- 
servative estimate. 
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There exists no suitable means of calculating the wing contribution 
to damping in the Mach number region from 0.80 to 1.2CL.and a smoothly 
faired curve through this region from subsonic to supersonic speeds is 
of no value except possibly for triangular wings having rounded airfoil 
sections and which are not in the presence of a wake from a canard 
surface. 

The means of calculating the isolated surface contribution to 
'damping by using references 5 to 8, 11, and 12 generally give values 
that are nonconservative but are in reasonable agreement with experiment 
if used within the limitations imposed by the theory. 

It should be noted that for computation of wing damping at super- 
sonic speeds, acceptable results are obtained only if rC^ and 

are both used since the effects are generally in opposite directions. 
Computations of the damping derivatives for swept wings at lower super- 
sonic Mach numbers by use of references 7 and 8 are extremely limited by 
the theory. 


CONCLUSIONS 


From a study of the pitch-damping derivatives measured in flight on 
about 33 models and full-scale airplanes of widely differing configura- 
tion and geometric characteristics, the following general conclusions 
are offered. 

1. No suitable means exist for calculating the wing damping in the 
transonic region and a smoothly faired line between the subsonic and 
supersonic values will, in most cases, not conform to the experimental 
d'ata and will probably give an unconservative result. 

• 

2. The damping contribution of unswept and swept wings of U5° or _ 
less sweep Is erratic in the transonic region and may be either positive 
or negative, leading to possible dynamic instability for tailless 

conf igurations . 

3. Configurations with triangular and swept wings having approxi- 
mately 60° sweep of the leading edge and rounded airfoil sections exhibit 
fairly smooth variations of damping in the transonic region. 

4. Calculations of the damping for configurations that have a hori- 
zontal tail surface appear to give the best agreement with the experi- 
mental data If the distance from the trailing edge of the mean aerodynamic 
chord of the wing to the tail Is used for computing the downwash lag. 
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5. Damping moments resulting from the down-wash on a rear lifting 
surface caused by the lift increments on a forward surface proportional 
to the pitching velocity and the rate of change of angle of attack may 
he appreciable if the forward surface is at a large distance from the 
center of gravity, as on a canard configuration. 

It is believed that no detailed conclusions are warranted because 
.of the unsystematic nature of the data and the effects of relatively 
unknown factors such as experimental accuracy, nonlinearities with lift 
coefficient and oscillation amplitude, and the effects of oscillation 
frequency. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va. 
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TABLE II 


REDUCED -FREQUENCY RANGE 
FOR DATA PRESENTED 


Model 

Range of 

ax:/2V 

1 

0. 0212 

0.0258 

2 

.0287 

.0455 

3 

.0166 

.0217 

it- 

.0118 

.0126 

5 

.0046 

.0053 

6 

.0119 

.0135 

7 

.0081 

.0108 

8 

.0050 

.0072 

9 

.0081 

.0097 

10 

.0122 

.0129 

ll 

.0161 

.0169 

12 

.0165 

.0178 

13 

.0084 

.0156 

14 

.0168 

.0207 

15 

.0164 

.0241 

l6 

.0151 

.0273 

17 

.0069 

.0134 

18 

.0089 

.0106 

19 

.0209 

-0295 

20 

.0103 

.0125 

21 

.0058 

.00 86 

22 

.0081 

.0100 


! 

Airplane 

Range of 

cdc/2V 

1 



2 

0.0449 

0.0539 

3 

.0155 

.0310 

4 

.0117 

.0162 
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Notes Coefficients are based on wing area and. wing chord of model 6. 
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Figure 5.- Geometric characteristics and damping coefficients 
wingless body and horizontal tail combination (model 5). 
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Figure 6.- Geometric characteristics and damping coefficients for an 
aircraft configuration having an unswept tapered wing with a 
hexagonal airfoil section (model 6 ). 
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Figure 7*- Geometric characteristics and damping coefficients for 
an aircraft configuration having an unswept wing with an 
NACA 65(gg) A0 °4.5 airfoil section (model 7)- 
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Figure 8.- Geometric characteristics and damping coefficients for an 
aircraft configuration having the X-l wing (model 8). 
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Figure 11.- Geometric characteristics and damping coefficients for an. 
aircraft configuration having a modified triangular wing (model 11) 
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Figure 13.- Geometric characteristics and damping coefficients for a 
supersonic aircraft configuration (model 13). 
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Figure 15.- Geometric characteristics and damping coefficients of an 
aircraft configuration having a 40° sweptback wing (model 15). 
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Figure 1 6 .- Geometric characteristics and damping coefficients of an 
aircraft having a 40° sweptback wing that incorporates inverse 
taper (model 1 6 ). 









MCA RM L52K20 


45 



Center-of-gravity location 
Model A 50. 8 percent H.A.C. 

Model B 60.8 percent H.A.C. 


Wing 


£, aq ft 
H.A.C . , in. 


2.31 

5.23 

24.10 


Horiaontal Tall 

4.00 


A 

5, aq ft 
M.A.C., in. 


2.19 

11.83 


All dimenalona are in Inchea. 



.6 .7 .8 .9 1.0 1.1 1.2 1.3 1.4 1.5 1.6 

K 


Figure l8.- Geometric characteristics and damping coefficients of a 
60 ° triangular wing missile configuration (model l8). 













NACA RM L52K20 


8.23 «*- 
-*|J- 1-2.04 


►12.89-H 

*\ n~ 3tS8 


Center-of-gravity location 
Model A Station 38.6 

Model B Station 40.4 


Wing 

AT 3.02 

ST »q ft 1.422 

i.A.C., in. 8.20 


All dimensions are in inches 




□ 



m 

IBMMMM 


m 


m 

MB— 


551 


—7 

Model B (low lirt) J 


/ 


Canard 

A 

8 , aq ft 


Model A (low lift) 

Calculated 

r — Experimental 


-Calculated — 
Experimental 


Increment due to canard C T and C T 5 ’ — 

0 J L q (high lift}') — 1111 

•7 .8 .9 1.0 1.1 1.2 1.3 1.4 1.5 1.6 

M 

Figure 21.- Geometric characteristics and damping coefficients for 
a missile configuration having essentially unswept wings and 
canard surfaces (model 22). 
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Figure 25.- Geometric characteristics and damping derivatives for the 

D-558-II airplane (airplane 4). 
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